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Abstract. We give a family of congruences for the binomial coefficients ('^_^^) 
in terms of multiple harmonic sums, a generalization of the harmonic numbers. 
Each congruence in this family (which depends on an additional parameter 
n) involves a linear combination of n multiple harmonic sums, and holds 
mod pSi+S . The coefficients in these congruences are integers depending on n 
and k, but independent of p. More generally, we construct a family of congru- 
ences J-2n.fc for {''^Zi) mod p^"+'^, whose members contain a variable number 
of terms, and show that in this family there is a unique "optimized" congruence 
involving the fewest terms. The special case k = 2 and n = recovers Wol- 
stenholme's theorem {'^^~^) = 1 mod p^, valid for all primes p > 5. We also 
characterize those triples {n,k,p) for which the optimized congruence holds 
modulo an extra power of p: they are precisely those with either p dividing 
the numerator of the Bernoulli number B„_2„_fe, or fc = 0, 1 mod p. 



1. Introduction 

In 1862 the Rev. J. Wolstenhome [25] noted the congruence that for all primes 

'2p - V 
p-l 

This result is now called Wolstenholme's theorem. Later it was found that the 
related congruence on harmonic numbers _ff„ := stating that for all primes 

Hp^i = (mod p2), 

which was discovered earlier (by E. Waring (55] in 1782 and again by C. Babbage 
[2] in 1819), is in fact equivalent to Wolstenholme's result. 

In the following 150 years, Wolstenholme's congruence has been generalized in 
many directions (see Mestrovic [18] for a survey). This paper considers generaliza- 
tions in two directions. The first direction treats a larger set of binomial coefficients, 
replacing 2p~l with kp—1. In 1900 Glaisher [5] showed that for all integers fc > 2, 



1 (modp^). 



(1.1) i^p-l)^^ 

holds for all p > 5. In 1999 Andrews [Ij extended Glaisher's theorem to g-binomial 
coefficients. 

The second direction obtains congruences modulo higher powers of p, by adding 
extra terms to the right hand side of Wolstenholme's congruence. In 2000 van 
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Hamme [lOj proved a result implying that for all primes p > 7, 

(1.2) Pr')-i+2/E7 



where :— X)J=i 7 ^'"^ harmonic numbers. Recently Mestrovic [17) showed 
that for any prime p > 11, 

(1.3) p^"^')^l-2p^i + V E - (-°d/) 

V p — 1 / 7 27 

' j=i ■' i<i<j<p-i 

This congruence involves the additional expression 

l<i<j<j9-l 

which is an example of a multiple harmonic sum, defined below. 

The main result of this paper is a simultaneous generalization and unification of 
these results, giving congruences for (^JT/) to arbitrary powers of p, which involve 
multiple harmonic sums. Our basic method formulates and exploits the existence of 
families of linear relations between certain multiple harmonic sums. The coefficients 
in our congruences are given by certain polynomials 6j.„(r), defined below. 

1.1. Main result. A composition is a finite ordered list {Ai,...,Aj} of positive 
integers. For ease of notation, we will denote by {Ai, . . . , A^}" the composition 

A = { Ai , . . . , Aj , Ai , . . . , Aj , . . . , Ai , . . . , Aj } 

consisting of a concatenated copies of {Ai, . . . , Aj}. 

Definition 1.1. For a composition A = {Ai, . . . , Aj}, and a positive integer n, we 
define the multiple harmonic sum 

1 



■Ai -Aj ■ 

n>ii>...>ij>l *1 ■ . ■ . ■ *j 

(By convention, if A = {Ai, . . . , Aj} and n < j, we set H{\] n) = 0) 

For k an integer, it has long been known (see e.g. [16 ) that the binomial co- 
efficient i^pZi") can be written as a linear combination of 'elementary symmetric' 
multiple harmonic sums -ff({l}-';p — 1): 

p-i 



''^^l^) =Y.^k-iypiH{{lV;p-l) 
P ^ j=0 



For a fixed non-negative integer n, we may truncate this equation after the 2n- 
th term and use estimates on the p-divisibility of the multiple harmonic sums 
H {{1}^ ; p — 1) to obtain the congruence 

(1.4) C^l!) ^ E(^ - l)Vi?({lF;P - 1) (mod p2«+3) 

which holds for all primes p >2n + 5. 

In Section 2 we show that the generating function for the elementary symmetric 
multiple harmonic sums H({l}^;p—1) (withp fixed) satisfies a functional equation. 
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which we use to derive identities involving these sums. These identities can be used 
to cancel some of the terms appearing in (|1.4|) . Our main result (like equation (|1.4|) ) 
gives for each non-negative integer n a congruence for the binomial coefficient {^pZi) 
mod p2ii+3^ involving multiple harmonic sums. However, our congruence involves 
only the first n elementary symmetric multiple harmonic sums (instead of the first 
2n such sums). 

The coefficients in our congruences are given by polynomials in k. We make the 
following definition. 

Definition 1.2. Let < j < rt be integers. We define the extremal polynomial 
bj.n{T) G Q[T] to be the unique polynomial of degree at most 2n+l satisfying the 
following conditions: 

(CI) fe,, „(r) = (T - 1)^" mod (T - l)"+i 

(C2) &j- „(r) = {-lyr^ mod t"+i 

A table of the extremal polynomials bj^n{T) for < j < n < 3 can be found in 
Section [L2l Now we can state our main result: 



Theorem 1.3. (Optimized Congruences) Letn > be a fixed integer. The extremal 
polynomials 6j.„(T) (j — 0,1, ... ,n) have integer coefficients, and the following 
hold: 

(1) For every prime p > 2n + 5 and every integer k >1: 

(1.5) (^p-l) ^ E h,Ak)P^H{{iy;p- 1) (mod ^^"+3). 

(2) If p = 2n + i is prime, then for every integer k >\, the above congruence 
holds (modp2"+2). 

(3) For every prime 3 < p < 2n + 1 and every integer fc > 1, the above 
congruence is equality: 

- j2hnik)p'H{{iy;p-l). 

Wolstenholme's congruence is the case n = 0, fc = 2 of Theorem 11.31 As another 
example, taking n — k — gives the congruence 

'3p - 1' 



p-\ 



= 1 +402piJ({l};p- 1) - 396?j2iJ({i}2;p_ 1) 
+ 216/iJ({l}^;p- 1) (mod/). 



Theorem 11.31 has four important features: 

(1) The coefficients bj^nik) appearing in the congruence (jl.Sp are independent 
of the prime p. 

(2) There are a large number of congruences for C'PZi) holding mod p , 
which involve the multiple harmonic sums H{{iy;p — 1) for 1 < j < 2n 
(see Theorem 13.31) . The congruences (|1.5p are optimized among these in 
only containing the terms iJ({l}^;p — 1) for 1 < j < n. 

(3) The restriction of the theorem to exclude certain small primes, depending 
on n, is necessary. The congruences may fail to hold (mod p2n+3-j j-pj. 
p — 2n + 3 (when 2n + 3 is prime), and also fail to hold for p = 2. 
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(4) The extremal polynomials hj,n{T) depend on n, and for fixed j their values 
at integers 6j,n(fc), which are the coefficients in the congruences, do not 
stabilize as n — > cxd (with the exception of bo^nik); see the tables in Section 
ll.2p . However they do satisfy many interesting congruences as n varies, 
which we address in Section |6l 
One may ask whether the coefficients bj n(k) appearing in the extremal con- 
gruences (|1.5|) are uniquely characterized by (|1.5p holding for all sufficiently large 
primes p; we conjecture this is the case, and discuss it in Section [T31 

An exceptional congruence will be a triple {k,n,p) such that the corresponding 
congruence given in Theorem ll.3l holds modulo an extra power of p. We characterize 
exceptional congruences for primes p > 2n + 3 as follows. 

Theorem 1.4. (Exceptional Congruences) Let n be a non-negative integer, p an 
odd prime. For all fc G Z, we have the following: 
(i) Suppose p > 2n + 5. The exceptional congruence 

'ftp - 1" 



p-l 



3=0 



b,^nik)p>H{{iy;p~ 1) (mod 



holds if and only if either fc = 0, 1 (mod p) or p divides the numerator of the 
Bernoulli number i?p-2n-3- 
(ii) Suppose p = 2n + 3. The exceptional congruence 

'kp - 1 



p-l 



^6,-„(fc)p'i?({l}^;p- 1) (modp2«+3) 

3=0 



holds if and only if k = 0,1 (mod p). 



We obtain Theorem 11.31 and Theorem 11.41 as special instances of a family TN,k 
of generalized Wolstenholme congruences, given in Theorem l3.3[ and with the J-N,k 
specified in Definition 13.41 The general congruence in the family J-N,k (which will 
hold for all sufficiently large primes p) is of the form 

(1.6) C^Ii) ^ E b,p'H{{iy;p~ 1) (mod 

where e G {1, 2} is chosen so that e = N (mod 2), and the coefficients bj are ratio- 
nal numbers. Each congruence in this general family is derived from (|1.4p . using 
linear combinations of identities among multiple harmonic sums (these identities 
are stated as Theorem l2.2[) . The optimized congruence (11.51) is distinguished as the 
unique congruence in the family J-2n,k satisfying — bn+2 = ■ ■ ■ = &2n = 0. 

1.2. The extremal polynomials 6j_„(r). In SectionlHlwe prove some interesting 
properties of the extremal polynomials 6j^„(T). Here we present data on these 
polynomials for small j, n in Table 1 below. 

Table 1 illustrates that &o,n(2^) — 1 for all n (this will be established in Section 
ini). While the definition of 6j^„(T) given earlier shows that it has degree at most 
2n + 1, in fact its degree is at most 2n (see Theorem 14. 7p . 

We next consider the coefficients bj^n{k) appearing in the extremal congruences 
given in Theorem 11.31 Values of the coefficients for fc = 2 and fc = 3 are given in 
Table 2 (resp. Table 3) below. 
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Table 1. Extremal Polynomials 6j^„(T) 








1 


2 


3 





1 








1 


1 


- r 






2 


1 


_y4 _^ - T 






3 


1 


2T« - 6T^ + 5r* 


-2T« + 6T^ ~ hT^ 


r« - 3T5 + 3T'' 






-T+ 1 


+T2 


-T3 



Table 2. Values 6j>(fc) for fc = 2 








1 


2 


3 


4 


5 





1 












1 


1 


2 










2 


1 


-2 


4 








3 


1 


14 


-12 


8 






4 


1 


-66 


68 


-40 


16 




5 


1 


382 


-380 


248 


-112 


32 



Table 2 shows that 61,2(2) — —2, 62,2(2) = 4, so that Theorem 11.31 reduces to 
Mestrovic's result ()1.3p in the case n = k = 2. 

Table 3. Values 6j,n(fc) for fc = 3 



n\j 





1 


2 


3 


4 


5 





1 












1 


1 


6 










2 


1 


-30 


36 








3 


1 


402 


-396 


216 






4 


1 


-6078 


6084 


-3672 


1296 




5 


1 


102786 


-102780 


66312 


-29808 


7776 



Tables 2 and 3 illustrate that for j > 1, the coefficients 6j,„(fc) which appear in 
the congruences (jl.Sp do not appear to stabilize as n — )■ cxd (holding j and k fixed) . 

1.3. Uniqueness conjecture. The statement of Theorem 1 1.31 raises an issue con- 
cerning whether the coefficients 6j,„(fc) above are uniquely determined by the con- 
dition that the congruences (jl.Sp hold for all sufficiently large primes p. Theorem 
14.11 asserts that there is a unique congruence of the form (|1.5p in the general family 
J'N,k- We believe that our family J-N,k actually contains all congruences of the 
shape (II. 6|) which hold for all sufficiently large primes p, in which case the individ- 
ual coefficients 6j.„(fc) above would be uniquely determined by (|1.5p . but we do not 
establish this. We formulate this belief as the following conjecture. 

Conjecture 1.5. (Uniqueness Conjecture) 

Let n > 0, k be integers, 60, 61, • • ■ , 6„, £ Q such that the congruence 

/fcp _ A ^ ^ h^p,Hi{iy;p- 1) (mod p^"+2) 
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holds for all sufficiently large primes p. Then bj = 6j,„(fc) for j = 0, 1, . . . , n. 

This conjecture might be difHcuh to resolve, in view of the fohowing consequence. 

Proposition 1.6. In the special case n — k ~ 1, the Uniqueness Conjecture \1.5\ 
implies that there are infinitely many primes p such that the numerator of the 
Bernoulli number Bp-3 is not divisible by p. 

This property of the BernouUi numbers is currently an open problem. 

Proof. Take n = k = 1, bo = bi — 1. These values do not agree with the values 
of the extremal coefhcients &o,i(l) ~ 1, ^i,i(l) = 0, so the Uniqueness Conjecture 
states that there are infinitely many primes p > 7 for which 

l=(P'])^l+pH{{l};p-l) (mod/). 



That is, there are infinitely many primes p for which p"^ \ H{{l};p— 1). Glaisher 
[7] showed that, for all primes p > 5, 

H{{l};p-l)^~^p' (modp^). 
The result now follows. □ 



In Section 3 we formulate a more general Strong Uniqueness Conjecture 13. 6[ 
which we show in Section 4 implies the Uniqueness Conjecture. 

1.4. Extensions of results. Theorem 11.41 connects certain exceptional congru- 
ences with p-divisibility of the numerators of certain Bernoulli numbers. Recall 
that theorems of Herbrand and Ribet say that for odd i in the range S < i < p ~ 2, 
p\Bp^i if and only if a particular piece of the class group of the cyclotomic field 
Q(Cp) has non-trivial p-part (see j23]. Sec. 6.3 for a precise statement). This raises 
the possiibility that our congruences may have an interpretation in terms of the 
arithmetic of cyclotomic fields. We do not currently know if that is the case, but 
in we will show that there is a parallel family of congruences related to p-adic 
L-function values at positive integers, which involves the 'power' multiple harmonic 
sums H{{j};p~l). In [20] we systematically investigate the structure of identities 
among multiple harmonic numbers that underly such congruences. 

1.5. Related results. The literature contains a vast collection of identities and 
congruences involving multiple harmonic sums and related sums, starting with work 
of Euler on harmonic numbers. Some of these involve binomial coefficient congru- 
ences (see Granville [9] for a survey). Our generalized congruences appear to have 
a structure not observed before, but for convenience we summarize some results 
from the literature for comparison. 

A number of congruences are known for the elementary symmetric multiple har- 
monic sums H{{lY;n) considered in this paper. In 1900 Glaisher ^ proved that 
for all odd r > 5 and all primes p > 7, 

Srip) := ?^H{{ir;p- 1) - H{{ir-^;p- 1) ^ (mod p% 

holds. In 1953 Carlitz [5] sharpened the congruence of Glaisher to show for all odd 
r > 5 and prime p > 7, 

ar\ - i{p-r){p-r-l){p-r-2) 5 
^-(^^=^ 24(p-r-3)(p-l)! 
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giving a relation with Bernoulli numbers. Along similar lines, Tauraso [21) shows 
that for any prime p > 7, 

H{{l};p-l) ^ 

There are also many congruences known that involve power sum multiple har- 
monic sums 

H{{r},p-l) = J2i- 

i=i 

Washington |24| provided a formula expressing these sums in terms of values of 
p-adic L-functions at positive integers. Congruences modulo p involving arbitrary 
multiple harmonic sums have been investigated by Hoffman |12| . 

Multiple harmonic sums appear in certain calculations in physics. Bliimlein ([3] 
discussed applications of multiple harmonic sums to quantum field theory. He com- 
puted a family of algebraic relations between multiple harmonic sums (A; n) which 
are independent of the upper limit of summation n. The Hopf algebra of quasi- 
symmetric functions has also been used by Hoffman to investigate both multiuple 
harmonic sums [13 and multiple zeta values [11 . 

Kontsevich [15^ considered the function 

n j 

U ' 

which is a 'twisted' version of a multiple harmonic sum. The expression Hn{x) 
may also be considered as a trancated version of the series — log(l — x). The 
corresponding sums for the dilogarithm J2^=i were considered by Elbaz- Vincent 
and Gangl [B], who determine functional equations satisfied by such series. These 
values are truncated forms of multiple zeta values, which is of great current interest 
(see Zagier [56 , Brown 4 ). 

2. Representing Binomial Coefficients in Terms of Multiple 

Harmonic Sums 

Our first object is to express binomial coefficients in terms of multiple harmonic 
sums. For n a positive integer, define the polynomial 

({n + l){T+\)-\\ 

^ / • j=i 

It will be useful to rewrite /„ in the form 

«T).n(i + !l±ir) 

Expanding the product above, we see that the coefficient of is (n -I- 1)^ times 
the j-th elementary symmetric function in 1, i, i, . . . , i. In other words, we have 

n 

(2.1) U{T)^Y.^n + iyH{{iy:n)T^ 

3=0 
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By convention, we take H{{l}^;n) — when j < or j > n. We also take 
H{(t>;n) = l. 

The polynomial /„ satisfies a functional equation relating T and — 1 — T. We 
can compute 

fni-l-T) = l.f[((n + l)i-l-T)+j) 

/ -|\ri " 



n'. 

i=l 



(In the third line we have made the substitution i — n + 1 — j). We expand this 
functional equation using (|2.ip to get 

Y,{n + iyH{{iy;n)T^ = (-1)" ^(n + - T)^" 

j>0 j>Q 



J2{-l)-+^n + iyH{{iy;n) ^ (^)t^ 

j>0 0<i<j ^ ^ 



i>0 Vj>' 



This holds identically in T. Equating the coefficient of on each side of the 
preceding equality and rearranging gives the following identity. 

Proposition 2.1. For all non-negative integers n, j , we have 

(2.2) {n + iyH{{iy;n) + 5](-l)"+'+i (') {n + iyH{{iy- n) = 

The sum above is finite (terms corresponding to i > n vanish). We thus 
have a family of linear equations (indexed by j) satisfied by the quantities {n + 
lYH{{iy-n),i = QX...,n. 

From the above we obtain a general set of identities expressing binomial coeffi- 
cients in terms of the H{yiy\ n). 

Proposition 2.2. Let n he a non-negative integer, k, Co,Ci,... indeterminates, 
and define 

(2.3) br-={k^iy + c, + {^ir+'+'j2(l 

Then the equation 

(2.4) (^(" + ^)-^)=f;6,(n+l)^i7({lF;r.), 
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holds identically in the indeterminates k, co,ci, . . .. Here the right side of ()2.4p is 
a finite sum, since i?({l}^; n) — for j > n. 

Proof. To begin, we use Equation (j2.ip to write 
^k{n + l) - 1^ 
n 



fn{k-l) 



= ^(fc-l)-'"(n + l)^i/({lH;n) 

Considering n fixed, we add to this equation a Unear combination of equations 
p.2p (where cj is the coefiicient of the equation indexed by j) to obtain the general 
formula. □ 



Remark 2.3. By making suitable choices of the parameters Ci in Proposition [2721 
we can arrange to have bj = for many j . Theorem 14.11 is obtained by optimizing 
this process. Special cases of the identities in this proposition were noted long ago. 
For example, Emma Lehmer f[16j. p. 360) used the identity 

'^"^^ =(-ifE(-i)Vi^({iF;A:). 

This particular identity is valid for all integers m, namely 

k 

= {-!)' J2i-iym^H{{iy;k). 

j=o 

However most identities produced above hold only for restricted values of m, namely 
m = k{n + 1) — 1, for fixed n. 

3. Congruences for {''pZi) Modulo Powers of p 

To obtain congruences for ('^^J^^), we will take k to be an integer and truncate 
the expansion of Proposition l2.2| with n — p — 1, p a, prime. To establish a bound 
on the error due to truncation, we need to establish some congruence properties of 
multiple harmonic sums. 

In the remainder of the paper, when p is understood to be a fixed prime, we 
sometimes abbreviate 

H{{iy) ■.= H{{iy;p-l). 

3.1. Congruence properties of multiple harmonic sums. Zhao ([22], Theo- 
rem 1.6) gives the following congruence involving multiple harmonic sums H ({ly ; p— 
1) and Bernoulli numbers: 

Proposition 3.1. Let p be a fixed odd prime, and j an integer with 1 < j < p — 3. 
Then we have 

i ^^f+i ' P 'ifj = mod 2 

(modp3),/,^l nrod2 

We provide additional congruence for H{{iy);p~ 1) for those j which are not 
covered by Proposition 13.11 

Proposition 3.2. Let p be a fixed odd prime, and j a positive integer. 
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fi) If i = P - we have H{{l}->\p — l) = \p (mod p^). 
(a) If j = p — 1, we have H{{iy);p — 1) = — 1 (mod p). 
(in) If j > p, we have H{{iy;p — 1) = 0. 

Proof. (i) We have 

~t 1 (p - 1) 
1 

1 v{p ~ 1) 
(p-1)! 2 

where in the last hne, we have used Wilson's theorem, which states that 
{p — 1)1 = —1 (mod p). 
(ii) We have 

1 



H{{1} 



p-i\ _ 



iP - 1)! 
= — 1 (mod p) 

(iii) In this case the defining sum is empty. 



□ 



3.2. A general family of congruences. We can now obtain our general family 
of congruences for (^.j/). The congruences are obtained from Proposition 12.21 bv 
truncation. We take some care to epress the error due to truncation in terms of 
Bernoulli numbers. 

Theorem 3.3 (General Wolstenholme-like Congruence). Let k be an integer. Let 
Co, Ci, . . . G Q &e given, and take bj 6 Q defined by 

(3.1) b, = {k-iy+c, + {-iy+'j2^^ 

i=0 

Fix an odd prime p which does not divide the denominator of any Ci, and let N be 
a non-negative integer. Define 

EM:=[^^'l)-Y.b,j^^H{{lY,p-l) 



P 

(i) If < N < p ~ A, we have 



En = ( ^—6^+1 + 6^+2 jp ^ (mod p ^ ) 

when N is even, and 



when N is odd. 
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(ii) If N — p — 3, we have 

bN+i 



E 



N 



bM+2 (mody^+^) 



(iii) If N = p — 2, we have 

Em ^ -bj,+,p^+' (mod 

(iv) If N > p ~ 1, we have 



En = 



In particular this implies that 



(3.2) 



kp — 1 
p-1 



N 



J2b,p'Hi{iy 

3=0 



mod p^'^^ if N < p — 4:, N even 
mod ifN<p-4:, N odd 

mod if N =p-'i 

modp^+i ifN=p-2 
mod if N >p-l 



(Congruence mod means equality) 

Proof. We apply Theorem 12.21 to obtain the equahty 



kp — I 
p-1 



Y,byH{{iy). 

Because n = p — 1 is even, the values of bi are independent of p. It follows that 



j=N+l 

We use Proposition 13. 2i statement (iii) to eliminate the terms in the above some 
with j > p, giving 

(3.3) En= b,p^Hi{iy) 

N+l<j<p-l 

The bj are p- integral (because we assume the Cj to be p-integral) . Propositions 13.11 
and 13.21 also tell us that the expressions II({1}^) are p-integral. We now separate 
into cases. 

Case (i-a): Suppose < iV < p — 5 and N is even. Due to the integrality of bj 
and i7({l}^), we have 

EM^bM+,p^+'H{{l}^+') + 6^.+2P^+2i/({l}^+2) 

+ 6Ar+3P^+3i?({l}^+^) (modp^+4) 

Proposition[S3]tellsusthati/({l}^+i) = ^(^jj^a)^ (mod p^) and H{{1}^+^) = 

-^^p (mod p). Also, Proposition 13. II or Proposition 13.21 give that iJ({l}^+'^) = 
(mod p). Combining these give 



E 



N 



N + 3 



N + 2. 



bN+i+bN+2 P''+' (mody^+4) 



as desired. 
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Case (i-b): Suppose < N < p — 5 and N is odd. Due to the integrality of bj 
and iJ({l}^), we have 

Er, ^ + (mod 

Proposition[33]tellsusthati/({l}^+i) = ^Bp_2-NP {mod p^) and H{{1}^+^) = 
(modp^). This gives 

En ^ (mod p^+% 

as desired. 

Case (i-c): Suppose N = p — A. Due to the integrahty of bj and H{{1}^), we 
have 

Ej, ^ 6^.+lp^+li^({l}~+l) + 6^+2p"+2h({1}^+2) (mod 

Proposition 13.11 tells us that i/({l}^+^) = -^^Bp^2-NP (mod p^). Proposition 
13.21 tells us that i/({l}^+^) = ^p (mod p^). Combining these gives 

Case (ii): Suppose N — p — In this case p.3p becomes 

i?JV = 5w+ip^+li7({l}P-2) + 6Ar+2p^+2i/({l}f-l) 

Proposition Ogives i7({l}P-2) = £ (^od p^) and i/({l}P-i) = -1 (mod p), 

so 

En ^ {^-^ - 6^+2) (mod 
Case (iii): Suppose N ~ p — 2. In this case (|3.3p becomes 

= 6Ar+ip^+li/({ir-l) 

Proposition [3?2] gives ^^({l}'"^^) = — 1 (mod p), so 

Em ^ -bM+iP^'^^ (modp^+2) 

Case (iv): Suppose N >p-l. Then shows that En = 0. 

□ 

Definition 3.4. We call the congruence p.2p the generalized Wolstenholme con- 
gruence associated with the data 

[fc, (co,ci, . . .),iV] , 

and we will say that bo, . . . ,bN are the generalized Wolstenholme coefficients asso- 
ciated with this data. 

We let FM.k denote the family of all generalized Wolstenholme congruences 
above, where TV, k are fixed and the other data varies. 
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Remark 3.5. For fixed fc, N, the family J-N,k has the structure of an affine hnear 
space over Q in the following way: ii B = {bo, . . . ,bN) and B' = (feQ,...,&'^) 
are the coefficients associated with the data [k, {co,ci, . . .), N], [k, {c'q, c[, . . .), N] 
respectively, and t G Q, then 



where the numbers on right hand side are the generalized Wolstenholme coefficients 
associated with the data 



In the next section we will focus exclusively on the case where = 2n is even. 
We will determine that the affine space of generalized Wolstenholme coefficients, 
for arbitrary k and N — 2n, has dimension n. The 'optimized' congruence will be 
uniquely determined among this n dimensional family. 

Let be a positive integer and set M = + 3 if is even, M = A^ + 2 if A^ is 
odd. It might be reasonable to expect that every congruence of the form 



that holds for all sufficiently large primes, actually comes from Theorem 13.31 We 
formulate this as the following conjecture: 

Conjecture 3.6. (Strong Uniqueness Conjecture) If k, m are integers with m > 0, 
and oo, . . . , a„ £ Q are such that 



holds for all but finitely many p, then this congruence arises from Theorem \3.S\ in 
the following sense: there are constants Co,ci,... 6 Q such that, if bQ,bi, . . . are 
defined by then we have a,; — bt for i — 0,1, ... , ip{m), where ip{m) = m — 2 

if m is even and ip{m) — m — 3 if m is odd (here, we take ai — for i > n). 

Proposition 14.21 below will show that this conjecture implies the Uniqueness 
Conjecture 11.51 

Now we consider some special cases of Theorem 13.31 In what follows we will 
write (co, . . . , for the sequence (cq, . . . , Cm, 0,0,...). 

As one example, fix a positive integer k and take the data [fc, ((fc — 1)^), 2] . This 
gives (6o, 6i, &2, ^3, ^4) — (1, k{k — 1), 0), so we get the congruence 

Corollary 3.7. For all integers fc and all primes p ^ 2,5, we have 



tB + {l- t)B' = {tbo + (1 - t%, ...,tbN + {l- t)b'j^). 



[fc, (teo + (1 - t)c'o, tci + (1 - t)c[, . . .), N] 




N 




ao + aipH{{iy;p- 1) + . . . + a„p"i/({l}";p - 1) (mod p™) 



This is a generalization of van Hamme's result (|1.2I) . 
Taking the data [2, (49, -18, 4), 6] gives 
{bo,bi,...,bG) = (1,14, -12,8,0,0,0), so we get the identity 




(mod p'^) 



Corollary 3.8. For all odd primes p, we have 




1 + UpH{{iy;p- 1) - 12p^H{{l}^;p- 1) 



+8p^H{{l}^;p-l) (mod/) 
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Corollaries 13.71 and 13.81 are special cases of Theorein l3.3[ corresponding to n = 1 
and n — 3, respectively. 



4. Optimized Wolstenholme Congruences 

We now state a version of our main result (Theorem II. 3|) . We show that when 
= 2n is even, it is always possible to choose the data (co, ci, . . .) so that bn+i = 
bn+2 = ■ . ■ = &2ji — 0. Moreover, this condition will uniquely determine the values 
of bj for < j < n. We will derive this result using Theorem 13.31 

Theorem 4.1 (Optimized Wolstenholme Congruences). Let integers k, n be given, 
with n > 0, and set N = 2n. Then there exist unique values &j_„(fc) G Q (j = 
0,1, ... ,n) with the following property: 

There exist co,Ci,... G Q .such that the generalized Wolstenholme coefficients 
bo,bi,...,bN associated with the data [k, {co,ci, . . .), N] satisfy bj — bj^nfk) for 
< j < n, and 

bn+l = bn+2 ~ . . . ~ &2n = 

Additionally, the Cj may be taken to be integers, so that the 6j.„(fc) are necessarily 
integers. 

In other words, for N = 2n, Theorem lS. 3\ produces a unique congruence of the 
form 

(4.1) P^"/) =bo + b,pHi{iy;p-l) + ... + bnP^H{{ir;p-l) (modp2«+3) 



p-1 

with bi £ Tj, which holds for all odd primes p ^ 2n + 
In the case p = 2n + 3 we have 

(4.2) i^^Zl) = bo + bipH{{lY;p-l) + . . . + b^p^H{{ir;p-l) (mod ^^"+2) 

Additionally, for odd p < 2n+ 1, we have 
'kp — 1 



p-l 



= bo + bipH{{iy;p- 1) + . . . + &„p"iJ({l}"; p-l) 



We will defer the proof of this result, first showing that it gives an implication 
between the two uniqueness conjectures we have made. 



Proposition 4.2. The Strong Uniquenss Conjecture \3.6\ implies the Uniqueness 
Conjecture \1.5[ 

Proof. Let n > 0, k be given, and suppose bo, . . . ,bn G Q are taken such that the 
congruence 



kp-1 
p-l 



l + ^bjp>H{{iy;p-l) (modp2"+2) 
i=i 



holds for all sufficiently large p. Suppose the Strong Uniqueness Conjecture 13.61 is 
true, and apply it with m = 2n + 2, to find that there are cq, ci, . . . G Q such that 



3 , 

3 



6, = (fc-l)^+c,+(-l)^+i^mc 



i=0 
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for J = 0, 1, ... , n, and 



i=0 

for j = 71 + 1, n + 2, . . . , 2n. Now, the uniqueness statement of Theorem 14.11 above 
says that bj — bj^nik) for j = 0, 1, . . . , n. □ 

For the proof of Theorem l4.1l we need some prehminary definitions and lemmas. 

Definition 4.3. Fix integers N, k, with N >0. Define VN^k C Z^+i to be the set 

j 

V / .7 \ 

C,; 



VN,k ■■= < {bo, ■.■,bN): 3co, ci, . . . G Z s.t. b, ^ {k ^ ly + Cj + (-1)'+^ ^ 

I 1=0 

In other words V^v.fe is the set of generalized Wolstenholme coefficients corresponding 
to integer data. We similarly define f. C Q^+^ to be 

V^^j, := I [bo, . . . , 6^) : 3co, ci, . . . e Q s.t. b, ^ {k - I)-'" + c, + (-1)^+1 ^ Qq| , 

the set of generalized Wolstenholme coefficients corresponding to rational data. 
The inclusion VN,k ^ ^jv fe induces an isomorphism 

VN.k ® Q = V^^^ 

of afhne spaces over Q. 

We have that VM,k is a coset of the subgroup 

Vn I (6o, . . . , ^at) : 3co, ci, . . . e Z s.t. 6, = + (-1)^+1 ^ (^^^ Q | C Z^ 

Note that Vn is independent of k. Wc then have the following: 
Proposition 4.4. For all integers N, k, with N >Q, we have VN,k = VN,i-k- 
This is not surprising, as 

^kp- 1\ _ /(I - k)p-V 
p-l)-[ p-1 

holds for all odd primes p and all k. 

Proof. As VN.k and V/v,i-fc are cosets of the same subgroup Vn < Z^+^, equality 
will follow if we can show V/v,fe H V/v,i-fc 7^ 

Taking co = ci = . . . = 0, we see (1, (-fc), (— fc)^, • . . , {-k)^) e Vat^i-a;. To see 
that this element is also in VN,k, set Cj = —{k — 1)^. Then 

J / ■ 



6, = (fc-l)^+c, + (-ir+i^M 







3 

'J 



{k-iy -{k-iy + {-iyY. -^y 

i=0 

i~ky 

□ 
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Lemma 4.5. For positive integers b, n, let Mn.b be the n x n matrix 

^b + i^ 



0<i,j<n 



Then det M„ 



Proof. Define n x n matrices i„ = ~ ((j-i))- unipotent lower- 

triangular and Un.b is unipotent upper-triangular, so both have determinant 1. We 
claim that A/„.t, — L„t/„.ti, so that detM„^b = det i„ det [/„,& = 1- Looking at the 
entry in slot (i, j), this equality reduces to the Vandermonde convolution identity 

b + i\ b 



□ 



Lemma 4.6. For all non-negative integers n, V2n is a free 'L-module of rank n. 
Additionally, the map n : V2n — >■ given by 

7r(6o, . . . , &2n) = {bn+l, ■ ■ ■ , &2n) 

is an isomorphism 

Proof. Here {bo, . . . , b2n) € Vzn is determined by the values of Cj for < j < 2n. 
We therefore have a surjective map : Z-^"+^ — V2n < taking (cq, . . . , C2„) 

to (ao, . . . , a2n) with 

In other words the matrix representing (p„ (with respect to the standard basis 
on 2^"+^) is given by 



\V / Q<i,j<2n 



If we identify row vectors of length 2n + 1 with the set of polynomials of degree 
at most 2n via the identification (oq, . . . , a2n) -H- ao + aiT + . . . + a2,iT^", then the 
j-th row of An is identified with the polynomial 

This means that the row span of A„ is contained in the set of polynomials f{T) sat- 
isfyint f{T) ~ —f{—l — T). Such polynomials can be written as Q- linear combina- 
tions of T-fi, (T+i)2, (T+i)2"-i. It follows that rank(l/2«) = rank(A„) < n. 

Next let i : Z" ^ Z2"+i, (xq a;„-i) ^ (xo, . . . , a;„_i, 0, 0, . . . , 0). We have 
Tro(pnoi{xo,...,Xn-i) = (yi, . . . ,2/„), where 



i=0 

Bv Lemma l4.51 this map is biiective. It follows that tt is surjective. Since rank(V2n) < 
n = rank(Z"), we must have that rank(V2ri) = n, and tt is bijective. □ 
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Proof of Theorem \4-l\ We need to show that there is a unique element of the form 

(5o,...,6„,0,0,...,0) eK,^ ,, 

and that 60, . . . , 6n G ^- It sufhces to show there is a unique element of this form in 
V2n,k- Because V2n,fc = (1, (fc— 1), (fc — 1)^") + V2n7 this is equivalent to showing 



there is a unique element of the form 

a = (ao, . . . , a„, -{k - l)"+\ 



,,-(fc-l 



\2n\ 



This is the same as finding an a g V^n with 7r(a) = ( — (fc — 1)"^^, . 
Such an a exists and is unique by Lemma 14.61 

That the values bj^n{k) agree with a polynomial in k will follow from Corollary 
431 below. ' □ 



^(fc- 1)^") 



We summarize the recipe for constructing the coefficients bj^n{k) given in The- 
orem 14.11 It will follow that these coefficients are interpolated by a polynomial 

Theorem 4.7. For fixed integers < j < n, the coefficients hj_n{k) given in 
Theorem \4.1\ are values of a polynomial 6j_„(T) at T = k, which is of degree at 
most 2n. This polynomial can be computed explicitly as follows. 
Let Mn be the n x n matrix 



M„= (-1) 
Let Dn be the (ri + 1) x n matrix 

Dn 



n-\-i 



I 



G)- 









a<i,j<n-l 



0<i<n 
0<j<n-l 



where S^ j is the Kronecker delta. Then M„ is invertible over the integers, and we 
have the matrix equation 

(4.3) 



/ &0,n(fc) \ 

bi,n{k) 




( {k- 
(fc- 


1)0 \ 


- Dn ■ ■ 


( (k-ir+' \ 

(fc- l)"+2 


\ bn,„{k) J 






1)" ) 







In particular this shows that bj^n{k) is given by a polynomial in k, of degree at most 
2n, having integer coefficients. 



Proof. The formula (|4.3p follows from the proof of Theorem 14.11 This formula 
implies that bj^k{k) is given by a polynomial in k, of degree at most 2n, having 
integer coefficients. □ 

Definition 4.8. For integers j, n, k, with j,n > 0, we let bj^n{k) denote the 
coefficients arising from Theorem 14.11 We call these extremal coefficients. We also 
denote by bj^n{T) the polynomial giving these coefficients, and call these extremal 
polynomials. By convention, we take 6j.„(T) — for n + 1 < j < 2n, and we say 
that bj,n{T) is not defined for j > 2n + 1. 

Theorem 14.71 provides a recipe for computing the extremal coefficients in terms 
of matrices involving binomial coefficients. In this sense they are like the Bernoulli 
numbers, which can be computed by a similar expression. 
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Proposition 14.41 savs that V2n.k — V2n,i-k- We may therefore make the substitu- 
tion k ■(r^ 1 — k in Theorem 14. 71 to get the foUowing: 

Corollary 4.9. Let k, n be given, and set N = 2n. Let Dn, Mn be as in the 
statement of Theorem \4.T\ Then 



(4.4) 



/ bo{k,n) ^ 
bi(k, n) 



( {-kf \ 

' i-ky ^ 



Dr, ■ m: 



V b,,{k,n) J V i-kr J V i-ky- ) 

We can combine Theorem 14.71 and Corohary 14.91 to obtain the foUowing charac- 
terization of the extremal polynomials 5j^„(T): 

Proposition 4.10. Fix integers j, n, with j < 2n. The extremal polynomial 
&j^„(T) e Z[T] is the unique polynomial of degree at most 2n satisfying the following 
conditions: 
ft) „(r) = (T - 1)^" mod (T - !)"+! 
fti) „(r) = {-Ty mod r"+i 

Proof. For n + 1 < j < 2n, 6j^„(r) — 0, and the result can be seen directly. For 
j < n, Corollary 14.71 shows that deg(6j^„(r)) < 2n, and that bj^niT) is equal to 
(T - ly plus a Z-linear combination of (T - . . . , (T - l)^".' This shows that 

ft), fti) are satisfied. Condition ftii) similarly follows from Corollary 14.91 

For uniqueness, the Chinese remainder theorem says that conditions fti), ftii) 
determine the residue class of bj^n{T) modulo 7^"+i(T — There will, a fortiori, 
be only one polynomial with rational coefficients and degree at most 2 n -I- 1 in this 
residue class. Our computations haye shows that this polynomial in fact has integer 
coefficients and degree at most 2n. □ 



Theorem 11.31 stated in the introduction, now follows from the combination of 
Theorem 14. 11 Theorem 14.71 and Proposition 14. 101 

5. Exceptional Congruences and Bernoulli Numbers 

We now investigate the situations under which the congruences (|4.1I) , (|4.2I) hold 
modulo some larger power of p than giyen by Theorem 14.11 We term these ex- 
ceptional congruences. In the case of Wolstenholme's theorem, we have that the 
exceptional congruence 

= 1 (mod p ^ 

,p-lj 

holds if and only if p divides the numerator of the Bernoulli number i?p_3 (this 
follows from the results of van Hamme [10] and Glaisher [7] ) . We establish a similar 
result, which shows that the congruences (|4.1|) . ()4.2|) hold modulo an extra power 
of p if and only if either p\Bp-2n-3, or p\k{k — 1). 

Let non-negative integers n, k be given, and choose co,ci, . . . G Z so that the 
generalized Wolstenholme congruence associated with the data [k, (cq, ci, . . .), 2n] 
is the optimal one, given by Theorem 14.11 (the Ci are not uniquely determined by 
this condition). Let &0; ^ii • • ■ i ^2n+2 be given by (|3.1|) . so that bj = foj.n(fc) for 
j = 0, 1, . . . , n, and bj — for j — n + l,n + 2, . . . , 2n. The values of &2n+i and 
&2n+3 will depend on the choice of the 
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Proposition 5.1. Define 

C„(fc) := {n + l)62ri+l + b2n+2 

Then, independent of the choice of Ci giving the optimal congruence, we have 

Cn{k) = fc"+l(fc - 1)"+^ 

Proof. First we show that the value of Cn{k) depends only on n and k, but not 
the choice of a. Let p : V2n+2,k V2n,k, (&o, • ■ • , &2n+2) ^ (&o,---,&2n) be the 
projection map (where V.^k is given by Definition 14. 3|) . From the construction of 
the spaces V2n+i.k, V2n,k, P is surjective. Lemma 14.61 savs that rank(V2n+2,fc) — 
n + 1, rank(V2Ti,fc) = n. If follows that U p^^{bo,bi, . . . ,bn,0, . . . ,0) is a Z 
torsor. Therefore, if we can exhibit (6q, . . . , &2„+2) (^Oi • ■ • : ^2n+3) G U such that 
{n + l)fo2n+i + ^2«+2 = {n + l)b2n+i + &2«+2, wc wih be done. 

If we take = Ci for i ^ 2n + 1, and C2„+i = C2n+i + 1, we will have that h\ = hi 
for i < 2n, &'2„+i = 62n+i + (^Y^)' ^nd 6i„+2 fo2n+2 - (^V^)- It follows directly 



that (n + 1)62,1+1 + ^2«+2 ^ {n + l)62n+i 

Next we show that Cn{k) agrees with a polynomial in k. Using the same process 
as in CoroUarv 14. 7[ we may solve for the data co,...,C2n to give the extremal 
congruence. We will use this data (with = for i > 2n + 1). We can then 
compute &2n+i, 62n+2 iu the following way: 
Let Mn be the n x n matrix 



Let A„ be the 2 x n matrix 



(-1) 



n-\-i 



(-1) 



i+i 



0<i j'<n-l 



2n+l<i<2n+2 
0<j<n-l 



where Sij is the Kronecker delta. Then M„ is invertible over the integers, and we 
have the matrix equation 



(5.1) 



02n+l 
b2n+2 



(fc 
(fc 



2^2n+2 



An ■ M- 



( {k 

(fc 



l)"+i 

iy-i+2 



\ 



J 



\ (fc-l)2" 

Moreover, 62n+i is equal to 
. , (fc- 1)^", so that 62„H 



2+1 IS 



This shows that 62,1+1, 62n+2 are polynomials in k. 
(k — l)2"+i plus a Z-linear combination of (fc — 1)""''^, . 
monic in k, of degree 2ri+l, and (fc — l)""''^|62n+i. Similarly, 62^+2 is monic of degree 
2n + 2, and (fc — 1)"^"^ |62n+2. It follows that C„(fc) = (n + l)62n+i +62,1+2 is monic 
of degree 2ri + 2, with (fc — l)^"+^|C„(fc). Moreover, C„(fc) is determined by the set 
V2n+2,fc, and Lemma lT^ savs that V2n+2,k — V2n+2,i-k- We may therefore make the 
substitution k ^ 1 — k to get the fc"+^|C„(fc). By the Chinese remainder theorem, 
fc"+^(fc — l)"+^|C„(fc). The only monic polynomial of degree 2n+2 which is divisible 
by fc"+i(fc - l)"+i is fc"+i(fc - l)"+i, so we conclude C„(fc) = k"+^{k - l)"+i. □ 

We now consider the possibility of extra powers of p in the congruences (|4.ip , 
(|4.2p . For all integers fc, n with n > 0, and all odd primes p, define 

fkp — 1 



E{k,n,p) 



\p-l 



Y,h,Ak)p'H{{lY) 
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Proposition 5.2. Suppose p > 2n + 5. Then 

2n + 3 

Proof. By Theorem 13. 3[ we have 

EE 62n+ip'"+'i?({l}'"+') + 52„+2p'"+'i?({l}'"+') (mod p2"+4) 
Usmg Proposition 13. 11 we can write 

Eik,n,p) EE -62„+iP^"+^-^^Bp-3-2„-W2P^"+^%^ 

2(2n + 3) 2n + 3 

= •^P-^-^" _|_ l)52„+i - 62«+2) 

n + 3 
n + 3 



Proposition 5.3. Suppose p = 2n + 3. Then 

E{k,n,p) EE -k"+\k - i)"+y"+2 (mod p2"+3) 

Proof. By Theorem 13.31 we have 

E{k,n,p) = Wi/"+'ff({l}'"+') + W2/"+'i?({l}'"+') 
Using Proposition 13.11 we can write 



□ 



El/- 7 N — &2n+l 2n+2 j, 2n+2 

E{k,n,p) = —^P -b2n+2P 

[p - l)b2n+l , \ 2„+2 
2 02n+l I P 

EE - {{n + l)b2n+l + b2n+2) 

EE -fc"+i(fc-l)"+V"+2 (modp2"+3) 



□ 



We can now state the precise conditions under which the congruences (|4.ip . 
(|4.2p hold modulo a larger power oip than is given by Theorem l4.11 The following 
Theorem is an immediate consequence of the preceding two propositions: 

Theorem 5.4. Let n > 0, k be integers, p an odd prime, and 5j^„(T) the extremal 
polynomials (characterized by Proposition \4.10^ . 

(i) Suppose p > 2n + 5. The congruence 

''^II) ^ E hn{kyH{{iy;p- 1) (mod p2«+4) 

holds if and only if either p|i?p_3__2n or k = 0,1 (mod p). 

(ii) Suppose p = 2n + 3. The congruence 

^^II) ^ E hn{k)p^H{{iy;p~ 1) (mod p^-+') 
holds if and only if k = 0,1 (mod p). 
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6. Properties of the Extremal Polynomials 

The extremal polynomials bj^n(T) satisfy many arithmetic and congruence rela- 
tions. 

Proposition 6.1. The extremal polynomials bj,n{T) satisfy the following proper- 
ties: 

(i) For all non-negative integers n, &o.n(T') — 1 
(ii) For all non-negative integers n, &„.„(T) ~ y''^{T — 1)" 
(Hi) For all non-negative integers j < 2n, {T — ly divides &j,„(T). 

Proof, (i) and (ii) follow by checking that the given polynomial satisfies the condi- 
tions of Proposition ITTOl (Hi) also follows immediately from Proposition 14. 101 □ 

For fixed j, the polynomials bj^n{T) depend on n. One exception to this is that 
bo,n{T) = 1 for all non-negative integers n. Examining the table in Section 1.2, we 
see that &i^„(r) does indeed depend on n. However, 6i,„(T) + b2,n{T) = — T for 
all n. 

This is the first in a family of equations giving linear combinations of the extremal 
polynomials 6j.„(T) (j varying) which are independent of n. 

Proposition 6.2. Let m be a non-negative integer. Suppose we are given f(T) — 
Ej=o '^f^^ ^ satisfying f{T) — /(— 1 — T). Then, for all non-negative integers 

n with 2n > m, we have 

m 

^a,6,,„(r) = /(T-l) 

Proof. Define 

m 

9{T) ^^ajbj^n{T) 
j=o 

By Proposition Oni we have 6j,„(T) = (T - 1)^ (mod (T - 1)"+^), so that 

m 

g{T) EE J2a,{T^iy 

j=o 

= /(T-1) (mod (T - l)"+i) 
Similarly, by Proposition OUl we have 6j- „(T) = {-Py (mod r"+i), so that 

m 

g{T) EE Y.^A-Ty 
j=o 
= fi-T) 

= f{P-l) (mod (T - l)"+i) 
By the Chinese remainder theorem, it follows that 

g{P) = f{P-l) (mod T"+i(r-l)"+i). Now, we just note that g(r) and /(T - 1) 
have degree at most 2n + 1, so they must be equal. □ 
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